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ABSTRACT 

We use high time cadence images acquired by the STEREO EUVI and COR 
instruments to study the evolution of coronal mass ejections (CMEs), from their 
initiation, through the impulsive acceleration to the propagation phase. For a 
set of 95 CMEs we derived detailed height, velocity and acceleration profiles and 
statistically analysed characteristic CME parameters: peak acceleration, peak 
velocity, acceleration duration, initiation height, height at peak velocity, height 
at peak acceleration and size of the CME source region. The CME peak accel- 
erations derived range from 20 to 6800 m s~ 2 and are inversely correlated to the 
acceleration duration and to the height at peak acceleration. 74% of the events 
reach their peak acceleration at heights below 0.5 R®. CMEs which originate 
from compact sources low in the corona are more impulsive and reach higher 
peak accelerations at smaller heights. These findings can be explained by the 
Lorentz force, which drives the CME accelerations and decreases with height and 
CME size. 

Subject headings: Sun: coronal mass ejections (CMEs), Statistic 
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Introduction 



Coronal mass ejections (CMEs) are sporadic 



ejections of ma 



the Sun with masses in the order of 10 1 3 -10 16 g (iVourlidas et al. 



>16 



in the range of ~100-3000 km s 1 (e.g. 



gnetized plasma from 



Yashiro et al 



2004 : 



201 0) and velocities 



Gopalswamy et al. 



20091 ). 



They are accelerated by magnetic forces in the solar corona and move outwards into the 
interplanetary space, where they may severely influence the space weather near the Earth. 



S everal case studies of the CME kinemat i cs fjGallagher et al.l l2003t IZhang et al.l 1200 ll . 



2004 



Maricic et al 



2004 : 



Temmer et al. 



2008, 



20101 ) showed that CMEs typ ically underg o 



three phases: a gradual evolution, fast acceleration and propagation phase ( IZhang et al. 
20011 ). In the gradual phase, the CME leading edge rises slowly with velocities of some 
ten km s _1 . At a certain height, the CME undergoes a strong acceleration. How fast and 
how long the acceleration takes place varies from event to event. After the main impulsive 
acceleration phase the CME propagates at almost constant velocity or shows a gradual 
acceleration/deceleration due to the interaction with the ambient solar w ind flow during 



the propagation in the interplanetary space (e.g. 



Gopalswamy et al. 



2000) 



Recent studies suggest that the main acceleration of impulsive CMEs occurs at low 
coronal heights, which are not observable in trad itional white light coronographic images 



( jGallagher et al. 



2003 



Temmer et al. 



2008 



20101 ). Thus, if we are interested in the origin 



of CMEs and their initial acceleration, we have to observe them from their initiation site 
close to the solar surface. To observe fast CME accelerations it is important to have 
image sequences with high time cadence during this important dynamical phase of the 
CMEs. Since fast and impulsive events severely contribute to our space weather, they are 
particularly relevant to study. 

Acceleration measurements of the earl y phase of the CME propagation for a 



statistical sample of 50 events were done by 



Zhang &: Derd (120061 ). who used coronagraphic 
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observations from 1.1 R Q to 30 R Q (LASCO CI, C2 and C3), before the LASCO CI 
coronagraph f ailed in 1998 due to the communication loss with the SOHO spacecraft. 



Vrsnak et al 



(120071 ) combined EUV images (SOHO EIT) with coronagraphic observations 
(MLSO Mark-IV K, LASCO C2 and C3) in order to track CMEs from their initiation site 
up to about 30 R Q . They analysed a sample of 22 events, which contain predominantly 
gradual CMEs. 



For our study we used Solar Terrestrial Relations Observatory (STEREO; 



Kaiser et al. 



20081 ) data, which provide high time cadence EUV imaging and coronagraphic observations 
up to 15 Rq with an overlapping field-of-view (FOV). Based on this data set we derived 
detailed CME kinematics and acceleration profiles for a sample of 95 impulsive CME events, 
which occurred during January 2007 and May 2010, representing the largest data set for the 
study of impulsive CME acceleration so far. In this paper statistics and correlation analysis 
of the kinematical and dynamical CME characteristics are presented. In our parameter 
study we focus on the CME peak velocity, peak acceleration, acceleration duration, height 
at peak velocity, height at peak acceleration, initiation height and the CME source region 
size. The relation of the CME characteristics to the associated flare, filament eruption, 
large scale EUV waves and magnetic topology will be presented in a forthcoming paper. 



2. Data and data reduction 



The Sun-Earth- Connection Coronal and Heliospheric Investigation (SECCHI) package 



([Howard et al.ll2008l ) onboard the twin STEREO spacecraft, STEREO-A ( Ahead) and 



Wuelser et al. 



STER EO-B (Behind), includes an Extreme Ultraviolet (EUV) Imager (EUVI; 
20041 ). two white light coronagraphs (COR1, COR2) and two white light heliospheric 
imagers (HI1, HI2), observing the Sun in different FOVs. We combined data from EUVI, 
COR1 and COR2, in order to study the kinematics of CMEs from their initiation close to 
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Fig. 1. — Time sequence of EUVI 171 A (first and second row), C0R1 (third row) and 
COR2 (bottom row) images of the CME observed on 5 April 2008 by STEREO-A. The 
curved line marks the identified CME leading edge, the straight line indicates the prop- 
agation direction. The CME kinematics derived for this event is plotted in Figs. and 
[HI All EUVI images are plotted in an x-range of [—1600", +1600"] and in an y-range of 
[-1470", +1730"], the selected COR1 range for both directions is [-3720", +3720"] and the 
COR2 range [—14 000", +14 000"]. Note that not all images available are shown. A movie 
for the whole event including all images available is included in the electronic supplement. 
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Fig. 2.— Same as Fig. □ but for the CME event observed on 23 May 2007 by STEREO-B. 
The CME kinematics derived for this event is plotted in Figs. [6] and |HJ The ranges in x- 
and y-direction are [-1580", +1580"] for EUVI 195 A, [-3720", +3720"] for COR1 and 
[-14 000", +14 000"] for COR2. 
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the solar surface up to a distance of about 15 R Q . EUVI observes the solar chromosphere 
and low corona in four different wavelengths in a FOV up to 1.7 R Q . To track a CME we 
mostly used images in the 171 A passband because of the high time cadence up to 75 sec, 
but in some cases it was only possible to track the CME in the 195 A observations, which 
on average have a lower time cadence (for most events 10 minutes, but in some cases as 
good as 2.5 minutes). The high time cadence enabled us to obtain detailed acceleration 
profiles especially at the beginning of the CME propagation. 

The two STEREO coronagraphs COR1 and COR2 observe the inner and outer solar 
corona in a FOV of 1.4 to 4 R Q and 2.5 to 15 R Q , respectively. The overlapping FOVs of 
all three instruments enabled us to connect the same structure in the observations by the 
different instruments. The time cadence of the COR1 observations is mainly 5 minutes but 
can be up to 20 minutes in some cases, the cadence of COR2 is 30 minutes. 

We started our study with a sample of 146 selected CME events. For our selection 
we preferred events that could be observed already in the low corona (i.e. by the EUVI 
instrument), in order to get insights in the early phase of the CME dynamics. In addition, 
we required that the source region of a CME could be identified on the visible solar 
hemisphere, for which we used the location of associated flares as an additional marker. Not 
each event fulfills both requirements (99 events could be measured in the EUVI FOV, for 
89 events a flare could be associated). In most cases the associated flares were weak. Five 
events were associated with an M class flare, 26 CMEs with a C class flare. The remaining 
flares were of GOES class B or below^ The final CME sample for the statistics presented in 
this paper comprises a set of 95 events, for which we could derive the full CME acceleration 
profile, i.e. the peak acceleration and the acceleration duration could be measured for each 

l We note that the CME events in this study occurred during the extreme solar minimum 
period (2007-2010), which explains the lack of associated high-energy flare events. 
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CME. 84 out of them were tracked starting from the EUVI FOV. 



All data were reduced by the SECCHI solar software routine secchi_prep. The 
EUVI image sequences were corrected f or differential rotat ion, and for weak events a 



normalizing-radial-graded filter (NRGF; 



Morgan et al. 



2006h was used. For the COR1 and 



COR2 observations a pre-event image was subtracted. In addition, a sigma filter and a 
normalization technique were applied to get better contrasts of the transient faint CME 
structures. For the measurements of the CME evolution, running difference and running 
ratio images were reconstructed. 



3. Methods and Analysis 

We derived the CME kinematics by measuring the position of the CME leading edge in 
EUVI, COR1 and COR2 images. In order to obtain reliable and reproducible measurements, 
an algorithm was developed which semi-automatically determines the position of the CME 
leading edge in EUVI, COR1 and COR2 running difference images. The leading edge of a 
CME appears as a bright front with a sharp intensity drop to regions outside the CME. 
This information is the basis for the algorithm developed to identify quasi-automatically 
the CME leading edge and its evolution in subsequent images. The algorithm works in 
the following way: The running difference images are contoured with brightness levels 
starting at very low, positive intensities. Then the distances of the pixels on the outermost 
contour with respect to Sun center are calculated. This process is repeated with contours 
at incrementally increasing brightness levels until the mean distance derived from two 
subsequent contour levels are separated by a sufficiently small, prescribed distance of each 
other. The position of the last determined contour is supposed to provide the location of 
the CME leading edge. The result can be manually corrected if necessary (e.g., if image 
artifacts are included in the contour, or if the CME front is not well defined). We note that 
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faint CMEs were mostly measured by visual identification of the CME leading edge, since 
here the algorithm fails. Figures [T] and |5] show examples of EUVI, COR1 and COR2 running 
difference image sequences of two well observed CMEs, where the identified CME leading 
edge and the determined propagation direction are indicated. We assumed a propagation 
along a straight line and defined for each event a direction, which crosses the outermost 
part of the CME front. The CME distance was derived by following the evolution of the 
identified CME leading edge along this main propagation direction. In a few cases a linear 
propagation could not be assumed due to deflection towards the heliospheric current sheet. 
For these events, the outermost part of the measured CME leading edge was manually 
selected in each frame. The distance was averaged over an angular extent of 10° of the 
determined leading edge and measured from the CME source region, which allows us to 
follow the CME propagation even if it deviates from radial direction. 

There are several factors affecting the determination of the CME leading edges and 
several possible sources for systematic errors. On the one hand, the CME itself can change 
its appearance with time, e.g., the front might become blurred which makes it difficult to 
track the same feature over several solar radii. On the other hand, we combine observations 
of different instruments with different angular resolution and the detection sensitivity varies 
over the FOV. In addition, stray light levels are different for COR1 and COR2, influencing 
the appearance of the observed white light feature. In order to estimate the average error 
included in our kinematical measurements, we analysed some selected CMEs in detail. For 
each of these events we followed four times the tracked CME leading edge in time in EUVI, 
COR1 and COR2 images. For each run, the scaling of the images was chosen differently, in 
order to account for the different visibilities of the outer (fainter) CME features. From the 
thus obtained height-time measurements, we derived the mean and standard deviation at 
each instant from the different measurement runs. The results from this procedure suggest 
an average error of 0.03 R Q for measurements in EUVI, 0.125 R Q in COR1, and 0.3 R Q in 
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COR2 data, respectively. 

Numerical differentiation of the height-time curve provides the CME velocity and 
acceleration profiles as a function of time. Due to the fact that errors in the height-time 
curve are enhanced by the first and second derivative, a smoothing and fitting method 
is used. We used a spline fitting procedure in which the measured CME height-time 
profile is subdivided into consecutive segments. All segments are then fitted by cubic 
splines and at their end points ('nodes'), the functions merge continuously and are twice 
continuously differentiable. The users's input to the fitting procedure is the number and 
position of nodes. The spline-fitted curve is then used as the basis to derive the velocity and 
acceleration evolution by subsequent numerical differentiation. In the course of fitting, we 
also estimated errors for the velocity and acceleration for each of the fitted segments. For 
each segment, the uncertainties on the polynomial spline coefficients are derived. The errors 
on the velocity and acceleration were then determined via Gaussian error propagation. 

The advantage of using this method is that it provides reasonable errors (especially 
for the acceleration) also in cases of a high time cadence, where the classical error 
estimates using two neighbouring points yield typically very high errors for the CME 
accelerations. Different sets of nodes giving similar fits for a certain height-time curve may 
re yeal considerably d ifferent velocity and acceleration profiles (for details see Appendix 



in 



Vrsnak et al 



20071 ). A fit was preferred over the other if the errors in velocity and 



acceleration were smaller, which basically confined the number of nodes to 5 to 7. 

To estimate the C ME source region size we used different methods, described in detail 



in 



Vrsnak et al. 



(120071 ): 



the distance between bipolar coronal dimming regions derived from EUVI 195 A or 
171 A observations; 
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• the footpoint separation of the associated eruptive filament/prominence observed in 
EUVI 304 A images; 

• the lengths of the flare ribbon brightenings observed in the chromosphere or transition 
region measured in EUVI 304, 171 or 195 A images. 

4. Results 

In this paper a statistical analysis of 95 impulsive CME events is presented. For 90 
events out of them a source region could be uniquely determined. Fig. [3] shows their position 
on the solar disc. Most of the events under study occurred close to the solar limb, and thus 
the influence of projection effects is small. Fig. [3] also reveals the transition from solar cycle 
no. 23 to no. 24, with the CME source regions changing from locations close to the equator 
at the end of cycle no. 23 to higher latitudes for CMEs already belonging to cycle no. 24. 
Fig. H] shows the distribution of the projected radial distance of the CME source regions 
to Sun center in units of R & . A clear peak at 0.9 to 1 R Q can be seen, which includes 49 
events (more than 51%). Only 21% have distances smaller than 0.7 Rq, and only one event 
occurred close to disc center. 

In Figs. [1] and [2] we show EUVI and COR image sequences of two sample CMEs from 
our data set. The two events occurred on 5 April 2008 and 23 May 2007, respectively, 
and could be well observed in all three instruments, EUVI, COR1 and COR2. In both 
cases 16 representative images are selected; the whole evolution of the events can be seen 
in the accompanying movies no. 1 and no. 2. In each running difference image the solar 
limb is overplotted and the identified CME leading edges are marked. We measured the 
CME propagation along a straight line (also plotted into the image) which originates at the 
CME source region. In both examples the observed CMEs change their shape during their 
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Fig. 3. — Positions of the CME source regions on the solar disc as observed from the ac- 
tual STEREO vantage points. Different colors and symbols represent the different years of 
occurrence during the minimum phase of solar cycle no. 23/24. 
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Fig. 4. — Distribution of the projected radial distances of the CME source regions from Sun 
center in units of R Q with a binsize of 0.1 R & . 
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Fig. 5. — CME kinematics for the events observed on 5 April 2008 (left) and 8 May 2007 
(right). The top panels show the height-time curve derived from EUVI (crosses), COR1 (tri- 
angles) and COR2 (asterisks) measurements together with the measurement errors. Note 
that these errors (0.03 R Q for EUVI, 0.125 R Q for COR1, and 0.3 R Q for COR2 measure- 
ments) may appear smaller than the plot symbols due to the large height range presented. 
The solid line represents the spline fit to the height-time curve. The middle and bottom pan- 
els show the CME velocity and acceleration profiles derived from numerical differentiation of 
the CME height time measurements of the spline fit (solid line). The velocity and accelera- 
tion values derived by direct numerical differentiation of the measurements points (symbols) 
as well as the error range derived from the spline fits (grey shaded area) are overplotted. 
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Fig. 6. — Same as Fig. [5] but for the events observed on 5 May 2010 (left) and 6 November 
2008 (right). 




Fig. 7. — Same as Fig. 0but for the events observed on 23 May 2007 and 13 February 2009. 
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Fig. 8. — Evolution of the CME velocity and acceleration in dependence of the CME height 
for the events occurring on 5 April 2008, 8 May 2007 (top), 5 May 2010, 6 November 2008 
(middle), 23 May 2007 and 13 February 2009 (bottom). Note that the height-time curves of 
these events are shown in Figs. [SHU The data points derived from the different instruments 
are marked as crosses (EUVI), triangles (COR1) and asterisks (COR2). The solid lines 
indicate the first and second derivative of the spline fits, respectively, which are surrounded 
by the estimated error range (grey area). 
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Fig. 9. — Distribution of the CME peak velocity (top) and peak acceleration (middle) derived 
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propagation until they reach the end of the COR2 FOV. Selected images from the different 
instruments, which observed the CME almost simultaneously, were overlaid to check if we 
really observe the same structure. This was possible for example for the EUVI image taken 
at 15:53:30 UT and the COR1 image taken at 15:55:00 UT in Fig. EJ 

The kinematical plots derived for the CMEs shown in Figs. [T] and [2] as well as for four 
more CME events from our sample can be seen in Figs. [5H71 The upper panels show the 
CME height, the middle panel the velocity and the bottom panel the acceleration evolution 
in time. Although all six events show the typical three phase kinematical behaviour 
(gradual initiation, acceleration and propagation phase), their kinematics reveal distinct 
differences. The event observed on 13 February 2009 (Fig. [7j) shows a very strong (a max 
~970 m s -2 ) and short acceleration phase of 9 min, whereas the events that occurred on 6 
November 2008 (Fig. [6]) and 23 May 2007 (Fig. [7j) exhibit long acceleration durations of 
50 and 72 min, respectively with peak accelerations of ~400 m s~ 2 and ~280 m s -2 . For 
the event observed on 5 May 2010 (Fig. [5]) an acceleration of 29 min and the smallest peak 
acceleration (a max ~250 m s~ 2 ) of the six events shown in Figs. [5HZI was measured. The 
remaining two events plotted in Fig. [6] (5 April 2008, 8 May 2007) reveal peak accelerations 
of 850 m s -2 and 690 m s -2 and acceleration durations of 26 min and 17 min, respectively. 

Fig. [8] shows for the six events plotted in Figs. [5H7] the CME velocity and acceleration 
profiles against height, respectively, together with the error ranges derived from the spline 
fits. Differences are remarkable, in particular at which height the peak acceleration takes 
place. The event observed on 6 November 2008 reaches its maximum acceleration at a 
height of 0.84 R & above the CME source region, whereas for the other five events the height 
at peak acceleration is distinctly lower. For instance the event observed on 23 May 2007 
was just at a height of 0.12 R & above the CME source region when it reached its peak 
acceleration of 280 m s -2 . 
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For each CME we derived several characteristic parameters: 

• peak velocity v max 

• peak acceleration a max 

• acceleration duration t acc 

• height at peak velocity h vmax 

• height at peak acceleration h amax 

• first height measured ho 

• source region size L 

The peak velocity v max and peak acceleration a max were derived from the spline fit 
of the velocity-time curve and the acceleration-time curve, respectively (see middle and 
bottom panels of Figs. EH?]). The measured velocity at the end of the COR2 FOV (~15 R Q ) 
could result from the combination of two different effects, the main and the residual 
acceleration. If the residual acceleration is positive, the CME velocity profile peaks later 
but if the residual acceleration is negative, the maximum velocity is reached at the end 
of the main acceleration phase. The residual acceleration can be positive for two reasons. 
On the one hand, in slow and gradual events the ambient solar wind flow causes a further 
acceleration. On the other hand it is also possible that th e CME accelerates continuously 



due to continuous energy release after the main phase (e.g. ICheng et al.ll2010l ). Since we are 
interested in the CME velocity corresponding to the main acceleration phase, we determined 
Vmax as that value of the CME velocity, when the CME acceleration has decreased to 10% 
of its peak value. 

The acceleration duration t acc =t acc _ en< i — t acc _ start was extracted from the acceleration 
profile, where t acc _ start and t acc _ end were definded as the times when the CME acceleration 
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profile is again at the 10% level of the peak value. The height at peak velocity h vmax and 
height at peak acceleration h amax were derived from the velocity and acceleration profiles 
against the height (see Fig. [SJ, h vmax being related to the 10% level of a max . The height h 
at which the CME was first observed can be understood as a rough estimate of the height 
of CME initiation, since we observed all CMEs from their origin in the low corona. 

The top panel of Fig. [9] shows the distribution of the CME peak velocities. We found 
a range for v max from 56 to 1279 km s _1 with a mean value of 526 km s _1 and a median 
of 460 km s" 1 . Th e mode of the distribution lies at 300-400 km s -1 . In comparison, 



Vrsnak et al. 



(120071 ) measured 22 CMEs, which occurred in a period between February 2002 
and January 2005, i.e. the maximum and decay phase of solar cycle no. 23, covering a range 
of 365 to 2775 km s _1 with a mean value of 940 km s _1 . These linear velocity values are 
higher than the CME peak velocities in our study, which covers only e vents that occurred 



during the extreme solar minimum. Indeed, 



Gopalswamy et al 



(120091 ). who derived the 



linear velocity for about 11000 LASCO CMEs that occurred between 1996 and 2006, found 
that the mean velocity values vary from 300 to 600 km s _1 during the solar cycle with an 
average value of 475 km s -1 for the whole sample. 

The distribution of the CME peak velocities v max derived in our stu dy is asymmetrical 



with a tail towards high velocity values. This coincides with the findings of 



Yurchyshyn et al. 



( 120051 ) for the distribution of the linear velocity of 4315 CMEs, which was fitted by a 
lognormal fit. The lognormal probability density function derived from our sample of v max is 
overplotted in Fig. (top panel) as a solid line. Mathematically, an independent variable x 
is lognormally distributed when its natural logarithm, ln(x), matches a normal distribution. 
A normal distribution is created by the sum of independent variables, whereas a lognormal 
distribution is created by the product of independent variables. In other words, if a variable 
is lognormally distributed, this hints at a multiplication of independent physical processes 
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underlying the distribution (ILimbert et al.ll2001t lYurchyshyn et al.ll2005l ). The probability 
function f(x) of the lognormal distribution can be written as: 



i 



2,71(7 X 



cxp 



(ln(x) — nY 
2^2 



(1) 



where \x is t he mean and a the s tandard deviation of the natural logarithm of x. With 



reference to 



Limbert et al 



( 120011 ) we term fj,* = e M and a* = e°~ as the median and 
the multiplicative standard deviation. Thus the confidence interval of 68.3% is given as 
[fi*/a*,fi* ■ a*}. From the lognormal fit to the v max distribution (Fig. EH top panel) we 
obtained ^=6.09 and o=0.47 corresponding to a median /x*=441 km s _1 with the bounds 
of the confidence interval at [276 km s _1 , 706 km s -1 ]. 

To check in which way projection effects influence our results we plot in Fig. [10] the 
CME peak velocity against the radial distance of the CME source region to Sun center. 
The correlation coefficient of these two parameters is c =0.011, i.e. very low and thus 
projection effects do not have a significant effect. This is probably due to the fact that there 
exists a distinct corre lation between the CME radial propagation and its lateral expansion 



( Schwenn et al. 



20051 ) . Thus, the effect of projection effect on the derived CME velocities is 



much less than it could be expected in the case of point sources. 

The middle plot of Fig. [9] shows the distribution of the CME peak acceleration a max for 
the whole sample. The a max values cover a wide range between 19 m s~ 2 and 6781 m s~ 2 . 
40% of the events have a peak acceleration >600 m s -2 , 19% show values >1000 m s~ 2 
and five events (observed on 30 May 2007, 3 Jun 2007 ~ 06:00 UT, 3 Jun 2007 ~ 09:00 
UT, 17 Aug 2007, 12 Feb 2010) reached a max values >2000 m s~ 2 . The distribution 
peak is well-defined at 200-300 m s~ 2 , the mean and median values are 756 m s -2 and 
414 m s~ 2 , respectively. For the distribution of the CME peak acceleration we again applied 
a lognormal fit, with the parameters /i=6.09 and o=1.04 corresponding to a confidence 
interval of [156 m s -2 , 1248 m s -2 ] around /i*=441 m s -2 . Since the majority of the a max 
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values are concentrated in the range between and 2000 m s 2 , we plotted the detail of the 
histogram for that range at the bottom of Fig. O 

The large range of a max valu es (19-6781 m s~ 2 ) s preading over two ord ers of magnitud e 



is similar to former studies from 



Vrsnak et al. 



Zhang fc Dere 



(120071 ). 40-7300 m s" 2 , and 
( 120061 ). 2.8-4464.9 m s~ 2 . The mean v alue of a max = 749 m s -2 in our study is comparable 



to the value from 



Vrsnak et al 



mean value in the sample of 



(|2007f). 840 ms 2 , but more than twice as high as the 
Zhang fc Derd ( 120061 ). 330 m s -2 . It is worth noting that on 



average the CME peak velocities derived from CMEs that occurred during the extreme 
solar mini mum (present study ) are considerably smaller than those during solar maximum 



studied by 



Vrsnak et al. 



( 120071 ) but the CME peak accelerations are similar in both samples. 
This may be an effect of the better time cadence of the STEREO instruments which enables 
us to reconstruct fast changes in the CME kinematics. 

The distribution of the acceleration phase duration t acc , which peaks at 20-30 min, is 
shown in Fig. [TT] together with the lognormal fit. More than 50% of the events have t acc 
values smaller than 30 min, however we measured also acceleration durations up to 8.6 
hours. Only one event (observed on 17 October 2008) has an acceleration duration >200 
min. Thus we show in the bottom panel of Fig. [11] a zoom into the range of 0-200 min. The 



smallest value we obtaine d is 4.5 min. 
than the values derived in Vrsnak et al 



le ari thmetic mean of 44.6 min is consider ably lower 



(120071 ) . 120 min, and 



Zhang fc Dere! ( 120061 ). 180 min 



which indicates that our sample contains m ainly impulsive event s. T he wide range o f t nrn o f 



three orders of magnitude was also found by 



Vrsnak et al 



(120071 ) and 



Zhang fc Dere! (120061 ). 



In both of these studies the distribution peaks at 0-50 min, consistent with our results. The 
distribution of t acc was fitted with the lognormal fit parameters /i=3.39 and cr=0.79 which 
corresponds to a confidence interval of [13.5 min, 65.4 min] around /z*=29.7 min. 



Fig. [12] shows the distributions for the various CME height parameters, ho, h t 
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and h amax . The distribution for ho, the height at which the CMEs were first observed, is 
shown in the top panel of Fig. [T2J This is a measure for the height above the solar surface 
at which the CME is initiated. However, we stress that h is not exactly the real CME 
initiation height but a rough measure for it, since it is affected by a) projection effects, 
and b) sensitivity issues, i.e. the measured ho is expected to be larger than the real CME 
initiation height in cases where the CME could not be identified from the very beginning. 
The measured ho distribution covers the range from 0.01 to 1.76 R & with a mean value of 
0.24 R Q and a median value of 0.14 R & . The maximum of the distribution is located at 
the low end at 0.1 to 0.2 R@. The large h values (> 0.6 Rq) are from events which could 
not be identified in the EUVI FOV but only in coronagraphic images. These may be due 
to events that really start from source heights > 0.6 Rq or may be related to observational 
restrictions in terms of sensitivity for faint CMEs. For our sample of 95 events this applies 
to 11 CMEs. 

The middle panel in Fig. [12] shows the distribution of the heights h vmax defined as the 
velocity values reached at the end of the CME main acceleration phase. The distribution 
starts at very low heights of 0.17 R & and extends up to 9.5 R & (i.e. close to the border 
of the COR2 FOV). 63% of the events are observed in the range 0-1 Rq. The mean and 
median values for h vmax are 1.46 Rq and 0.78 Rq, respectively. 

The bottom panel in Fig. [12] shows the distribution of the heights h amax at which the 
CME accelerations reach their maximum. The mean of the h amax distribution is 0.53 Rq, 
the median 0.26 R & . The peak in the h amax distribution lies between 0.2 and 0.3 Rq, 
and 74% of the events have h amax smaller than 0.5 Rq. This means that most of the 
CMEs under study reach the acceleration at very low heights above the solar surface and 
emphasizes the importance of CME observations in the low corona in order to study the 
main acceleration phase. 
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Table 1 gives an overview of all the statistical CME parameters derived: we list the 
minimum and maximum value, the arithmetic mean with the standard deviation, the 
median together with the mean absolute deviation (mad) and the two fit parameters of the 
lognormal fit (fi and a) for v max , a maX} t acc , h , h amax and h vmax . 

In order to identify general characteristics and relationships, which intrinsically 
describe the evolution of CME eruptions, we correlated the various CME parameters 
derived. The correlation plots are shown in Figs. [131 US! an d HHJ All the correlations are 
plotted and calculated in logarithmic space. We found several CME parameters which 
revealed a significant correlation with the CME peak acceleration a max . The top panel in 
Fig. [13] shows the scatter plot of a max against t acc revealing a high negative correlation with 
a correlation coefficient of c = —0.84, i.e. CMEs with larger peak accelerations have shorter 
acceleration durations. This fi nding fits well with the results obtained in former studies 



covering smaller CME samples ( iZhang fc Dere 



200 



Vrsnak et al. 



of the regression line in Fig. [13] is in accordance with the results of 
who found a slope of —1.14. 
the following power law relat 



2007) . The slop e (-1. 09) 



Vrsnak et al. 



(120071 ) 



who found a slope of —1.14. The dependence between a max and t acc can be described with 

ionjfl 



a = 1Q 4 -23 4-1.09 (2) 
u>max J-U ^acc • 

The height h amax at which the CME reaches its maximum acceleration, and the peak 
acceleration a max are also anti-correlated, with c = —0.54 (Fig. [T3J middle), i.e. CMEs 
which are accelerated at lower heights reach higher peak accelerations. Between a max and 



2 A11 relations between the different CME parameters are calculated in the same units as 
used in the scatter plots (Figs. [131 [151 and IT6]) . i.e. m s~ 2 for a max , km s -1 for v max , min for 
tacc sud Rq for ho, h vmax and h amax . 
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mimimum 


maximum 


arithmetic mean ± 
standard deviation 


median ± mad 


/i ± a 


v max [km s" 1 ] 


56 


1279 


526 ± 263 


460 ± 160 


6.09 ± 0.47 


a max [m s~ 2 ] 


19 


6781 


757 ± 1034 


414 ± 246 


6.09 ± 1.04 


t acc [min] 


4.5 


516 


44.6 ± 60.4 


29.0 ± 14.5 


3.39 ± 0.79 


ho [Rq] 


0.01 


1.76 


0.24 ± 0.29 


0.14 ± 0.08 




hvmax [Rq\ 


0.17 


9.5 


1.56 ± 1.82 


0.78 ± 0.42 




hamax \R&\ 


0.04 


2.90 


0.53 ± 0.64 


0.26 ± 0.12 





Table 1: We list the statistical CME parameters derived: the peak velocity v max , the peak 
acceleration a max , the acceleration phase duration t acc , the height h , where the CME leading 
edge could be identified for the first time, the height at peak velocity h vmax , the height at 
peak acceleration h amax . Minimum value, maximum value, arithmetic mean with standard 
deviation, median with the mean absolute deviation (mad) are derived from the whole data 
set of 95 events, /i and a are derived from the lognormal fit to the distribution. 
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hamax we found the following dependence: 

a - in 2 - 32 h-° m ft) 

Umax — lu ' L amax- \°) 

There is also a distinct correlation between v max and a max , c = 0.58 (Fig. [THl bottom), 
which is not surprising, since CMEs which have stronger accelerations are also capable 
of reaching higher peak velocities. Between these two parameters we found the following 
power law dependence: 

a = in- - 60 y 1 - 21 (A) 

Umax iu u max- X*) 

In contrast to the peak accelerations, the CME peak velocities do not show a significant 
correlation with t acc , h vma x and h ama x (plots not shown). 

Fig. IT41 shows as a test the CME peak velocity v max against t a ccO ma x/2, which 
corresponds to modeling the CME acceleration by a triangular profile. We obtain a high 
correlation coefficient of c = 0.89. For comparison, the regression line (solid line) is plotted 
together with the 1:1 correspondence (dashed line), which are nearly parallel but shifted 
against each other by ~0-100 km s -1 . This difference can be attributed to a residual 
acceleration of the CME not captured by the simple triangular profile assumed. For the 
linear regression line we found the following relation: 

Vmax = 0.94 t acc ^p + 102.1, (5) 

with the same units (km s -1 ) used on both sides of the equation. 

The acceleration duration t acc against h vma x and h am ax also shows a power law 
dependence with correlation coefficients of c=0.76 and c=0.58, respectively (Fig. 115)) . Their 
power law dependence can be written as: 

t = 10 151 h om (6) 

t = 10 L73 h - 52 (7) 

L acc 1U 11 amax- V I 
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This means that CMEs which accelerate over a longer time also reach their maximum 
acceleration and velocity at larger heights. 

In Fig. [16] we show correlations of the height ho, at which the CMEs were first observed 
against h amax , t acc and a max . A strong correlation between ho and h amax is found with a 
correlation coefficient of c=0.82, i.e. CMEs which start at lower heights also reach their 
peak acceleration at lower heights. Between ho and t acc a weak correlation of c=0.43 
is found, i.e. CMEs which originate at low heights tend to accelerate more impulsively. 
Between ho and a max an anti-correlation with c = —0.46 was found, i.e. CMEs starting 
at lower heights in the solar corona reach larger peak accelerations. The relations can be 
expressed as: 

ho = 10 -0 - 40 h° a ^ ax , (8) 
h = lO" 1 - 56 (9) 
h = lO ' 25 a-° a f . (10) 

For 78 out of the 95 events it was possible to measure at least one feature to estimate 
the CME source region size. Because there are general differences between the different 
measurement methods we considered each feature separately. The flare ribbon length of the 
associated flares did not show any distinct correlation with the CME peak acceleration and 
velocity but the footpoint distances of the erupting filaments and the size of the coronal 
dimming do. Fig. [17] shows the filament footpoint distances (measured for 24 events) 
together with the linear extent of the associated coronal dimmings (measured for 10 events) 
against a max , h and t acc . The negative correlation between the source region size L and 
a-max (c = —0.50) and the positive correlation between L and t acc (c =0.38) indicates again 
that CMEs which originate from compact sources reach high er peak accele r ation s and have 
shorter acceleration durations consistent with the findings of IVrsnak et al.l (120071 ). Between 
the CME source region size L and initiation height ho a positive correlation was found with 
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the correlation coefficient of c =0.47. 

5. Summary and Conclusions 

In the following we summarize the most important findings of our study based on a 
sample of 95 impulsive CMEs observed in STEREO EUVI, COR1 and COR2. 

1. The histograms of the CME peak velocity v max , the CME peak acceleration a max 
and the CME acceleration duration t acc can be approximated with a lognormal 
distribution. 

2. We found a wide range of values for a max (19-6781 m s~ 2 ) and t acc (4.5 min-8.6 h). 

3. Most of the events (74%) reach their peak acceleration at heights <0.5 R e . 

4. a max and t acc are negatively correlated, c = —0.84. 

5. a max and the height at peak acceleration h amax are negatively correlated, c = —0.54. 

6. t acc and the height at peak velocity h vmax are positively correlated, c =0.76. 

7. t acc and h amax are positively correlated, c =0.58. 

8. t acc and the CME initiation height h are positively correlated, c =0.43. 

9. h and h amax are positively correlated, c =0.82. 

10. h and a max are negatively correlated, c = —0.46. 

11. The CME source region size L and a max are negatively correlated, c = —0.50. 

12. L and h are positively correlated, c =0.47. 

13. L and t acc are positively correlated, c =0.38. 
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Based on the assumption that the Lorentz force is the main driver of the eruption, it 
can be assumed that magnetic energy is transformed into kinetic energy and thus 

,2 B 2 



pv 
~2~ 



< 



2/i 



(11) 



or 



v < 



MoP 



1/2 



VA, 



(12) 



with B the magnetic field strength within the CME body, /xq the magnetic permeability 
in vacuum, va the Alfven velocity, p the plasma density and v the CME velocity. Eq, 
states that the CME velocit y cannot be larger than the Alfven velocity in the erupting 



structure (e.g., 



Vrsnak 



2006|). 



The Lorentz force density can be expressed as: 



/ = 3 x B 



1 

Mo 



(B ■ V).B — V 



B 



2\ 1 



(13) 



with j the current density. Inserting the approximation V « 1/L, we obtain the following 



order of magnitude estimate (see also 



Vrsnak et al. 



20071): 



a < 



B 2 



(14) 



2p pL 2L 

with L the characteristic length scale over which the magnetic field varies, which can be 
approximated by the CME source region size. Eq. [TH shows that the acceleration is not 
only governed by the Alfven velocity but is also dependent on the size of the erupting 
structure. Initially compact CMEs (small L, large va) will reach higher accelerations. 
These considerations coincide with our observational findings (summary item 10 and 11), 
i.e. the inverse proportionality found between a max and L (Fig. [TTJ, top) and also a max and 
h (Fig. dni bottom), which can be used as an alternative estimate of the source region size. 
We tested the relation between the CME peak acceleration, the size of the erupting flux 
rope and the Alfven velocity implied by Eq. [TH in that we derived for each event the ratio 
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v maxl '(2L) and v^ nax /(2ho) (with v max < va, see Eq. H2"]) . and correlated these quantities 
with a max . The resulting correlation coefficients lie in the range of 0.6-0.7 (Fig. [TSJ), i.e. 
are higher than the correlations of the individual CME parameters, which supports our 
interpretation. 

The distinct anticorrelation found between ho and a max with c ~ —0.5 is in addition 
also related to the stronger magnetic fields in the lower corona, which in turn are related to 
larger Lorentz forces providing the driving force for the CME acceleration. We note that 
we interpret the parameter ho, defined as the height where a CME is first observed, to be a 
measure of the initiation height of the erupting CME (and thus, to some extent, also to be 
a measure of the original size of the erupting structure). This is of course somewhat critical, 
since in gradual and faint CMEs we expect that we cannot really observe a CME from its 
very inititiation site (due to the limited sensitivity) but only further out when sufficient 
mass was accumulated at the CME front to be observable. Indeed, in 11 out of 95 events 
we could not find signatures of the erupting CME in the EUVI but only in the COR1 FOV. 
This may be a real effect, i.e. some CMEs in our sample started at heights > 0.7 R®, but 
may be partly also biased by the fact that the CME actually started at lower heights but 
was to faint to be observed. However, even if we underestimated the CME initiation height 
by using the height ho where the CME was first observed, the distinct correlations that we 
obtained between h and other characteristic CME quantities such as a max , v max and t acc , 
in line with the interpretations in terms of the Lorentz force elaborated above, support it 
to be a very useful quantity for CME initiation studies. 

Since during its propagation the size of a CME increases with height, we expect from 
Eq. [14] also that the acceleration decreases with height. This is consistent with our findings 
of a distinct anticorrelation between a max and h amax (summary item no. 5, see also bottom 
panel of Fig. [T3|) . The distinct inverse proportionality that we derived between a max and 
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tacc (summary item no. 4, see top panel of Fig. [T31 explains the relative small range of 
CME velocities (1 order of magnitude) despite the large ranges, over which a max (2 orders 
of magnitude) and t acc (3 orders of magnitude) vary. 

Taking into account that the acceleration time cannot be shorter than the Alfven wave 
signal transit time, we get the order-of-magnitude relation: 

tacc > — 15 

va 

i.e., CMEs originating from compact sources accelerate more impulsively, consistent with 
our summary item no. 13 (Fig. [T71 bottom) and no. 8 (Fig. [161 top), where we interpret ho 
as another measure of the source region size. 

The distributions of the CME peak velocity, CME peak acceleration and acceleration 
duration show a lognormal behavior (summary item no.l). Such distributions are created 
by the product of several independent variables. We suggest that the energy of the CME 
is dependent on at least two major variables. On the one hand, it depends on the amount 
of initially stored magnetic energy, which is available for transformation into other forms of 
energy. Furthermore, the CME energy depends also on the "transmission coefficient", i.e., 
the percentage of the initially stored energy which will be transformed into kinetic energy 
of the CME. 
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Fig. 10. — CME peak velocities v max against the projected radial distance to Sun center in 
units of Rq. 
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Fig. 11. — Distribution of CME acceleration duration t acc together with a lognormal fit. The 
bottom plot shows a zoom-in restricted to t acc < 200 min. 



-34 - 



40 



i2 30 

c 
> 

° 20 
<u 

e 

c 10 



o 

0.0 0.5 1.0 1.5 2.0 
80 1 I 




30 



25 

CO 

5 20 
> 
<y 

o 15 
i_ 
cu 
r> 

E 10 

=> 
C 

5 


0.0 0.5 1.0 1.5 2.0 2.5 3.0 

^OmOX ["sun] 



Fig. 12. — Distribution of h , i.e. the height at which the CME front could be first identified 
(top), height at peak velocity h vmax (middle) and height at peak acceleration h amax (bottom). 
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Fig. 13. — CME peak acceleration against acceleration duration (top), height at peak ac- 
celeration (middle) and peak velocity (bottom). The solid lines are linear regression lines 
to the data points with k the slope and d the ^/-intercept; c gives the correlation coefficient. 
Note that a max as well as t acc , h amax and v max are plotted on a logarithmic scale and that 
the fits and correlations were also determined in logarithmic space. The same holds for Figs, 
and I 
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Fig. 14. — CME peak velocity against the acceleration duration multiplied by half of peak 
acceleration, together with the regression line (solid). For comparison also the 1:1 corre- 
spondence is plotted (dashed line). 
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Fig. 15. — CME acceleration duration against height at peak velocity (top) and height at 
peak acceleration (bottom). 
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Fig. 16. — Height h at which the CME front was first detected against the height at peak 
acceleration (top), acceleration duration (middle) and peak acceleration (bottom). 
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Fig. 17. — Estimates of the CME source region size against the CME peak acceleration 
(top), the height of the first CME observation (middle) and acceleration duration (bottom). 
Crosses indicate the distance of the associated filament footpoints, diamonds the spatial 
extent of the coronal dimming region. 
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Fig. 18. — Top: Relation of Eq. [14] with the CME source region size L (filament footpoint 
distance, dimming extent). Bottom: Same relation but with the height ho at which the CME 
was first observed instead of L. 



-41 - 
REFERENCES 

Cheng, X., Zhang, J., Ding, M. D., & Poomvises, W. 2010, ApJ, 712, 752 
Gallagher, P. T., Lawrence, G. R., & Dennis, B. R. 2003, ApJ, 588, L53 
Gopalswamy, N., Lara, A., Lepping, R. P., et al. 2000, Geophys. Res. Lett., 27, 145 
Gopalswamy, N., Yashiro, S., Michalek, G., et al. 2009, Earth Moon and Planets, 104, 295 
Howard, R. A., Moses, J. D., Vourlidas, A., et al. 2008, Space Sci. Rev., 136, 67 
Kaiser, M. L., Kucera, T. A., Davila, J. M., et al. 2008, Space Sci. Rev., 136, 5 
Limbert, E., Stahel, A. W., & Abbt, M. 2001, Bioscience, 51, 341 
Maricic, D., Vrsnak, B., Stanger, A. L., & Veronig, A. 2004, Sol. Phys., 225, 337 
Morgan, H., Habbal, S. R., & Woo, R. 2006, Sol. Phys., 236, 263 

Schwenn, R., dal Lago, A., Huttunen, E., & Gonzalez, W. D. 2005, Annales Geophysicae, 
23, 1033 

Temmer, M., Veronig, A. M., Kontar, E. P., Krucker, S., & Vrsnak, B. 2010, ApJ, 712, 1410 

Temmer, M., Veronig, A. M., Vrsnak, B., et al. 2008, ApJ, 673, L95 

Vourlidas, A., Howard, R. A., Esfandiari, E., et al. 2010, ApJ, 722, 1522 

Vrsnak, B. 2006, Advances in Space Research, 38, 431 

Vrsnak, B., Maricic, D., Stanger, A. L., et al. 2007, Sol. Phys., 241, 85 

Wuelser, J., Lemen, J. R., Tarbell, T. D., et al. 2004, in Society of Photo-Optical 

Instrumentation Engineers (SPIE) Conference Series, Vol. 5171, Society of Photo- 
Optical Instrumentation Engineers (SPIE) Conference Series, ed. S. Fineschi & 
M. A. Gummin, 111-122 



Yashiro, S., Gopalswamy, N., Michalek, G., et al. 2004, Journal of Geophysical Research 
(Space Physics), 109, 7105 

Yurchyshyn, V., Yashiro, S., Abramenko, V., Wang, H., & Gopalswamy, N. 2005, ApJ, 619, 
599 

Zhang, J. & Dere, K. P. 2006, ApJ, 649, 1100 

Zhang, J., Dere, K. P., Howard, R. A., Kundu, M. R., & White, S. M. 2001, ApJ, 559, 452 
Zhang, J., Dere, K. P., Howard, R. A., & Vourlidas, A. 2004, ApJ, 604, 420 



This manuscript was prepared with the A AS IATjrpC macros v5.2. 



